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Effective hadron models commonly require the computation of functional determinants. 
Q ' In the static case these are one-loop vacuum polarization energies, known as Casimir 
Q\ ■ energies. In this talk I will present general methods to efficiently compute renormalized 
one-loop vacuum polarization energies and energy densities and apply these methods to 
construct soliton solutions within a variational approach. This calculational method is 
\J2 ■ particularly useful to study singular limits that emerge in the discussion of the classical 
Casimir problem which is usually posed as the response of a fluctuating quantum field to 
| externally imposed boundary conditions. 
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Introduction 



Many models in hadron physics originate from integrating out the more fundamental 
degrees of freedom like, for example, quarks. This change of field variables requires 
i efficient tools to compute functional determinants like 

^ ! Det (ifi+TiQi-m) , (1) 

in the Nambu-Jona-Lasinio model [|IJ for quark flavor dynamics 0. Here <J>j denote back- 
ground fields that couple to the internal symmetries via the generators IV Functional 
determinants are highly ultraviolet divergent and thus need to be regularized, and even- 
tually renormalized by perturbative counterterms. This is particularly elaborate when 
the perturbative expansions become invalid, as for example for soliton configurations 0. 
These configurations are usually static, §i(x) = &i(x ), such that the determinant becomes 
proportional to the vacuum polarization energy 

IEK-^ 0) )- (2) 

n 

Here u n are the eigenvalues of the single particle Hamiltonian in the presence of <&i(if) 
while o^ ** denote the eigenvalues in the case that the $j assume their vacuum values. 

1 Talk presented at the 2 nd intl. workshop Effective Theories of Low Energy QCD. 
2 Heisenberg-Fellow 
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In this talk I will describe tools to unambiguously regularize and renormalize such vac- 
uum polarization energies and energy densities starting from the energy density operator 
in quantum field theory. 

This presentation is based on work with E. Farhi, N. Graham, R. L. Jaffe, V. Khemani, 
M. Quandt, and M. Scandurra. The publications @, ||, ||, |7], |j of these collaborations 
should be consulted for further details. 



Method 

In quantum field theories energies and energy densities are computed as renormalized 
matrix elements of the energy density operator T 00 . Here I will present the method to 
unambiguously compute such a matrix element when the fluctuating quantum field is cou- 
pled to a classical background. The method is based on expressing this matrix element 
in terms of a Green's function with appropriate boundary conditions. Then the energy 
density is given by a sum over bound states plus an integral over the continuum scatter- 
ing states. Ample use will be made of analytic properties of scattering data especially to 
deform momentum integrals along a cut on the positive imaginary axis. To regulate the 
ultraviolet divergences of the theory, which corresponds to eliminating the contribution 
associated with the semi-circle at infinite complex momenta, the leading Born approxima- 
tions to the Green's function are subtracted and later exactly added back in as Feynman 
diagrams. These diagrams are then regularized and renormalized in ordinary Feynman 
perturbation theory. 



Formalism 



Consider a static, spherically symmetric background potential a = a{r) with r = \x\ in n 
spatial dimensions. The symmetric energy density operator for a real scalar field coupled 
to a is 



+ 



(-V 2 + 



m + air 



+ IV 2 (0 2 ) 



(3) 



with the spatial derivative term rearranged for later use of the Schrodinger equation to 
evaluate the expression in brackets. The "vacuum" is the state of lowest energy in 
the background a. The "trivial vacuum" is the state |0) of lowest energy when a = 0. 
The vacuum energy density is the renormalized expectation value of Too with respect to 
the vacuum \ fl), {fl\T 00 (x)\fl) ien , which includes the matrix elements of the counterterms. 
The energy density only depends on the radial coordinate r since a spherically symmetric, 



e(r) 



2tt"/ 2 



r-^TooOr)!^ 



(4) 



The wavefunctions factorize in radial functions <f>i(t,r) and angle dependent spherical 
harmonics. The Fock decomposition for the radial functions reads 



n-l 

r 2 



dk 



—y 



7TUJ 
1 



r 2 



if>t(k, r) e~ iwt a t (k) + ^(jfe, r) e iujt a\(k) 
^■(r)e- iw '4 fj +^(r) e^a^ 



(5) 
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This decomposition contains scattering states with uj =\Jk 2 + m 2 and bound states with 
uj£j = ^m 2 — njj. The total angular momentum assumes integer values £ = 0, 1, 2, ... in 
all dimensions except for n — 1, where I = and 1 only, corresponding to the symmetric 
and antisymmetric channels respectively. The radial wavefunctions ip are solutions to the 
Schrddinger-like equation 

- V" + £ (y - \) (v + |) V> + - fc V = ( 6 ) 

where i/ = £— 1-Mf • In each angular momentum channel, the wavefunctions are normalized 
to satisfy the completeness relation 
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oo 

dkif>i(k,r)if)t(k,r') +^^ej(r)^ ej (r') = 6(r - r') . (7) 



Then the standard equal time commutation relations for the quantum field <fi yield canon- 
ical commutation relations for the creation and annihilation operators, [ae(k), a\,(k')] = 
S(k — k')5u' and [atj, = 5jj'5u>- All other commutators vanish. The vacuum \Q) is 
annihilated by all of the a^{k) and a^. The matrix element (|j) can now be computed by 
inserting eq. (0) into eq. (0): 



e(r) = X> / — u^(M^M+Et^ 

r /"oo ^7j^ 1 

+J A- E ^ / — # (*, r)^(fc, r) + £ — ^ (r) 



e (0) (r) + e CT (r). 



Here N e is the degeneracy factor, D r = f — ^^J, ecT( r ) is the counterterm contri- 
bution, and e^(r) indicates the subtraction of the energy density in the trivial vacuum. 

The scattering state contribution is identified with the Green's function by defining 
the local spectral density 

p e (k,r) = -G t (r,r,k), (9) 



where 



G<(r,r, dg (fc + , e)2 _ g2 ~E fc2 + ^ . ( 10 ) 



o 



so that for real k 

= Im {kG e {r, r, k)} = Re{ Pi (k,r)} . (11) 

The ie prescription has been chosen so that this Green's function is meromorphic in the 
upper half-plane, with simple poles at the imaginary momenta k = in^j corresponding to 
bound states. For real k, the imaginary part of the Green's function at r = r' is an odd 
function of k, while the real part is even. Hence eq. (§) can be expressed as a contour 
integral in the upper half-plane. The contribution from the semi-circular contour at large 
\k\ with \m(k) > must be eliminated. Subtracting sufficiently many terms in the Born 
series from the Green's function yields a convergent integral. Then I add back exactly 
what I subtracted in the form of Feynman diagrams. The integrand (§) has a branch cut 
along the imaginary axis, k e [im, -Moo], and simple poles at the bound state momenta 
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k = iKj. The corresponding residues cancel the explicit bound state contributions. The 
discontinuity along the cut is hence all what is left to be considered 



e r 



m 7T 



4(t 2 - m 2 



-D r 



N 



[Pe(it,r)} N + Y / 4U r ) + e CT 



A/ 



= e r 



+ E 



-FD 



i=l 



where 



\pt(k,r)] 



pe{k,r) - pf'(k,r) - pt'{k,r) ... - p^'(k,r) 



(12) 



A? 



G*(r, r, k) - Gf (r, r, k) - G\"> (r, r, k) 



Gf } (r,r,A;)l . (13) 



The superscript (j) indicates the term of order j in the Born expansion. Subtraction 
of the free Green's function Gf\r,r,k) corresponds to subtracting e^(r) above. The 



potentially divergent pieces are precisely identified 



with Feynman diagrams e 



(0 
FD V 



which are regularized and renormalized using standard methods. When combined with 
the contribution from the counterterms ecT(^) they yield finite contributions to the energy 
density (for smooth backgrounds). 



The Radial Green's Function 

A variety of solutions to eq. (pi) is distinguished by different boundary conditions: 



free Jost solution wi(kr) 

Jost solution fe(k,r) 

regular solution </>i(k,r) 

physical scattering solution ipi(k,r) 



(-l^ypr [J u (kr) + iY v (kr)\ 
ij m MM = 1 



The physical scattering solution is normalized with respect to the Jost function, F^{k) 
that is obtained as the ratio of the interacting and free Jost solutions at r = 0, 



F e (k) = lim 



fi(k,r) 



wg(kr) 



(14) 



In particular, two regular solutions emerge: (pi has a simple boundary condition at r = 0, 
so that it is analytic in the upper half fc-plane; ipi has a physical boundary condition at 
r — > oo, corresponding to incoming and outgoing spherical waves. 
The Green's function has the simple representation 



(15) 



where r> (r<) denotes the larger (smaller) of the two arguments r and r'. The poles of 
Ge(r,r', k) occur at the zeros of the Jost function, which are the imaginary bound state 
momenta. These are the only poles of eq. (|T0|) in the upper half-plane and, since the 
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two functions in eq. ([H]) and eq. (|I5D obey the same inhomogeneous differential equation, 
they are indeed identical. 

Although Gg is analytic in the upper half-plane, fg and 4>g contain pieces that oscillate 
for real k and exponentially increase or decrease when k has an imaginary part. Actually, 
only the case r = r' is interesting, cf. eq. fll3D . Then the product fgcpg is well-behaved. 
This motivates to factorize the dangerous exponential components^, 



f e (k, r) = w e (kr)gi(k, r) and <pe{k,r) 



(_*)-"+* ht(k,r) 



2u we(kr) 

where wg is the free Jost solution introduced above. With these definitions, 

h £ (k,r)g e (k,r) 



(16) 



G e (r, r, k) 



2vg t {k, 0) 



(17) 



The definition of hg does not just remove the free part. Instead, it enforces the cancellation 
of W£ in the Green's function. After analytically continuing to k — it, g £ {it,r) obeys 

cf&it, r) = 2t&(tr)g' £ (it, r) + v(r)gt(it, r) (18) 

with the boundary conditions lim r _ ! . 00 <%(it, r) = 1 and lim,.-,^ g'^it, r) = 0. A prime 
indicates a derivative with respect to the radial coordinate r. Using these boundary 
conditions, one integrates the differential equation numerically for gi(it,r), starting at 
r = oo and proceeding to r = 0. Similarly, hi(it,r) obeys 



h'l(it,r) = -2t&{tr)ti £ {it,r) 



a{r) - 2t 



2 d &( T ) 



dr 



T=tr- 



h e (it,r) . 



(19) 



The factors in eq. fllBD were chosen to yield simple boundary conditions: hi(it, 0) = and 
h' e (it, 0) = 1. The numerical integration for h starts at r = and runs to r = oo. For real 

r, 

^( T ) = ~]n[wt(iT)] (20) 

is real with lim r ^ 00 ^(r) = 1, so the two functions hi(it,r) and gi(it,r) are manifestly 
real. They are also holomorphic in the upper half fc-plane and, most importantly, they 
are bounded according to \gi(k,r)\ < const, and \he(k,r)\ < const. [w/(i+]fc|r)]. Thus the 
representation of the partial wave Green's function in terms of gt and hg is smooth and 
numerically tractable on the positive imaginary axis. 

The computation of the Born series, eq. (|13D , is also straightforward in this formalism. 



The solutions to the differential equations eq. ([18]) and eq. (0) are expanded about the 
free solutions, 



gi{it,r) 
he(it,r) 



l + g P(it,r)+gf\it,r) + ... 
2vrl v {tr)K v {tr) + h[ l \it,r) + hf\it,r) + 



(21) 



where the superscript labels the order of the background potential a. The higher order 
components obey inhomogeneous linear differential equations such that a is the source 



3 For n = 1 and n = 2, the case of t = is somewhat different 



5 



term for ag^ is the source term for g( 2 \ and so on. Substituting these solutions in 
the expansion of eq. ( |T7D with respect to the order of the background potential finally 
yields the Born series for the local spectral density 



\pt(H,r)] 



N 



he(it,r)g e (it,r) 
2ug t (it,0) 



(22) 



J N 



Thus I have available a computationally robust representation for the Born subtracted 
energy density e(r), cf. eq. (|T2|). 



Feynman Diagram Contribution 

To one-loop order, the Feynman diagrams of interest are generated by expanding 

-l" 



(0|T 00 (x)|0)~|Tr 



T 

j_ ... 



-d 2 



m 



a 



(23) 



to order N in the background a. Here T x is the coordinate space operator corresponding 
to the insertion of the energy density defined by eq. @ at the spacetime point x, and the 
trace includes space-time integration. The Feynman diagrams are obtained in ordinary 
perturbation theory, thus the matrix element in eq. (p3| ) is evaluated between the trivial 
vacuum state, which is annihilated by the plane wave annihilation operators. The energy 
density operator has pieces of order a and a 1 : T x = T^ ' + T^ 1 '. The computation of 
its vacuum matrix element is most conveniently performed in momentum space. The 
relevant matrix elements are 



(k'\f^\k) 



i(k'—k)x 



k°'k° + k' -k + m 2 } and (k'\f^\k) = a( 



x)e 



i(k'—k)x 



(24) 



Here I will explicitly consider the contributions to (0|T o(a;)|0) that are linear in a. The 
first contribution of this order comes directly from T^> 



;Tr 



-d 2 



m 



2 x 



\a{x) 



d d k 



(2vr) d k 2 - m 2 ' 



(25) 



This local contribution is ultraviolet divergent for d = n + 1 > 2. It is canceled identically 
by the counterterm in the no-tadpole renormalization scheme. An additional contribution 
at order cr originates from Tj°) and the first-order expansion of the propagator, 



:Tr 



f(0)_ 



(26) 



-d 2 — m? x —d 2 — 

At 0(a) the renormalized Feynman diagram contribution to the energy density becomes 



4^(r) + 6 C T(r) 



C d r n ~ l 



d: 



d-l 



(2tt; 



C(l - C)q 



m 2 + ai-()q- 



2-d/2 



(27) 



with Cd = 27r^ i r(2-|)/[r(^)(47r) d / 2 ] and a(q) = 2Tr5(q°)a(q) is the Fourier transform of 
the (time independent) background field. This piece is finite for d = n < 4 and does 
not contribute to the total energy because it vanishes when integrated over space. The 
extension to higher order Feynman diagrams is straightforward. 
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Total Energy 

The total energy is simply the integrated energy density of eq. (|12| 



E[a] 



e{r)dr 



(28) 



Both, the t integral and the sum over channels in eq. (0), are absolutely convergent. This 
is a consequence of deforming the momentum integral in the upper-half planeQ. Thus the 
order of integration can be interchanged 



E 



cr 



00 dt 



7T 



poo N 

Jo i=i 



(29) 



where the total derivative term has integrated to zero. As already explained, a sufficient 
number, N, of Born approximations to the local spectral density pg(it,r) must be sub- 
tracted to render the t integral convergent. These subtractions are then added back in 
as the contribution to the total energy from the Feynman diagrams. Combined with the 
contribution from the counterterms this gives a finite result, 



N poo 

£4d + £ct=/ dr 



N 



£ e Fb(r)+e CT (r) 



Li=l 



(30) 



In practice, E[a] is more efficiently computed directly from the perturbation series of 
the total energy. First, I'd like to recall that 



roc d d 

2 1 dr \p(it,r)] N = - \hiF t (it)] N = - [\ng e (it,0)] N 
is validQ for Re(t) > 0, cf. Appendix A of Ref. 0. This allows me to write 

r°° dt t 



d 



(31) 



2tt \/t 2 — m' 



mt,o)] N +j2^i+E CT . 



The real function (3f (t,r) = In gf (it, r) is determined by the differential equation 



ffl(t, r) - Iftit, r)f + 2t&(tr)P' e (k, r) + a(r) 







(32) 



(33) 



with the boundary conditions lim^oo (3(t, r) = lim^oo (3'(t, r) = 0. In [/3e(t, 0)} N , the first 
iV Born terms must again be subtracted. They are obtained by iterating the differential 
equation ([33|) according to the expansion of Pi(t,r) in powers of a. 
To make contact with previous work j|, ^ [7], |TU|] , 



E\a\ = Y,N t 



N 



J2 E 



w 

FD 



E 



CT 



(34) 



4 For real momenta, k = it, this amounts to performing the momentum integral before the radial 
integral. 

5 In general, the case Re(t) = lm(fc) = causes uncontrollable oscillations at large r: for real k the 
integral dr [p(k, r)] N does not exist and the integrated local spectral density cannot be related to the 
phase shift. 
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observe that along the real axis the phase of the Jost function is the scattering phase 
shift, 

ilnF e (k)=iln\F e (k)\ + 5 e (k). (35) 

Equations (j32|) and (j34|) are proven identical by first noticing that for real k \Fe(k)\ and 
Si(k) are respectively even and odd functions and then computing the momentum integrals 
along the branch cut k G [im, +ioo]. 



Soliton Formation in a D=l+1 Chiral Model 

The described method is well suited to efficiently and unambiguously compute vacuum 
polarization energies in renormalizable quantum field theories. Then the total energy, 
i.e. the sum of the classical and vacuum polarization energies, is a functional of the 
background field. Varying this background field maps an energy surface. The existence of 
a local minimum on that surface indicates the existence of an energetically stable solution 
to the equation of motion, a solitonF]. Of particular interest are models that do not 
contain soliton solutions at the classical level such that solitons get stabilized by quantum 
corrections. 

Now I would like to consider this idea in the framework of a simple chiral model in 
.0 = 1 + 10. The realistic D = 3 + 1 case is more difficult and a discussion is presented in 
Ref. 0. In this two-dimensional model a two-component boson field = (0i, 02) couples 
chirally to a fermion \1/ that come in Nf (equivalent) modes: 

C = - ■ dy + *i W - G (0i + *750 2 ) } • (36) 

^ 71=1 



where the potential for the boson field 

A 



v{4>) 



2av 212 



A 



av 3 (0i — v) + const. (37) 



contains a term (proportional to a) that breaks the chiral symmetry explicitly in order 
to avoid problems stemming from (unphysical) infra-red singularities that occur when 
the vacuum configuration would be determined via the naive treatment of spontaneous 



symmetry breaking JTT|. In this manner it is guaranteed that the VEV is given by (0) = 
(v, 0). Here the counterterm Lagrangian is not presented explicitly. It is determined such 
that the quantum corrections lead to a vanishing tadpole diagram for the boson field. 
Considering only the classical contribution does not support a stable soliton soliton. 

In the limit that the number of fermion modes becomes large with v 2 /Nf ~ 0(1) only 
the classical and one fermion loop pieces contribute. In the following I will only consider 
that limit, i.e. E tot = E c \ + Ep. The fermion contribution can be split into two pieces 
Ep = O vac + E va \. The valence part E va \ is given in terms of the bound state energies 
such as to saturate the total fermion number that is fixed to be iVV. The vacuum piece 
is computed according to the formalism described in the preceding section: 

1 b - s ' r°° dk d / \ 

£v40] = -- E (kl - Gv) - j o — (u k - Gv) - (5 F (k) - 5^(kj) , (38) 

6 The minimum itself is not necessarily a soliton because the space of variational parameters is limited 
and the exact soliton might have even lower energy. 
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Figure 1: The maximal binding energy as a function of the model parameters as obtained from the 
Ansatz (|40|) in units of Gv; a = a/G 2 and A = X/G 2 . 

which is obtained from Eq ( |51D by employing Levinson's theorem. Note the overall "-" 
sign for fermions and recall that the single particle spectrum is not charge conjugation 
invariant. Furthermore, Sf denotes the sum of the eigenphase shiftsQ. The subtraction 

DC 12 „oo 

S^(k) = y o dx (y 2 - <f 2 {x)) (39) 

that renders E vac finite contains both first and second order Born approximants in the 
fluctuations of <p about ((f)). The first order is unambiguously fixed by the no-tadpole 
renormalization condition and the second order by the chiral symmetry. 

Having established the energy functional I now consider variational Ansdtze for the 
background field that turn this functional in a function of the variational parameters. As 
an example I assume 

fa + ifo = v {1 - R + Rexp [i-K (1 + tanh(Gvx/w))}} (40) 

that introduces width (w) and amplitude (R) parameters. For prescribed model param- 
eters (G,v,etc.) the energy must be minimized with respect to w and R. The resulting 
binding energy B = E tot — Gv is shown in Fig. fj. Even though the Ansatz (|40"D may 
not be the final answer to the minimalization problem, B is definitely negative. Thus a 
solitonic configuration is energetically favored showing that indeed quantum fluctuations 
can create a soliton that is not stable at the classical level. 



Casimir energies 

As noted earlier, the presented method to compute vacuum polarization energies is not 
limited to the case of smooth background fields. It is particularly interesting to employ 
singular background fields to imitate boundary conditions of the fluctuating quantum 
field. It is known for some time [0 that the vacuum polarization energy diverges when 



the fluctuating field is constrained by boundary conditions. Ad hoc schemes for their 
removal have been proposed However, the method of renormalization in continuum 



quantum field theory represents the only physical way to treat these divergences. 



7 The eigenchannels are labeled by parity and the sign of the single particle eigenenergies. 
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Dirichlet Points in One Space Dimension 

The simplest example to be considered is that of a massive scalar field <p(t, x) in one 
dimension, constrained to vanish at x — —a and a. The standard approach, in which the 



boundary conditions are imposed a priori, gives an energy [14 

E 2 {a) ■■ 



m 

T 



2a 

7T 



77T 



dt^m — r 
e iat - 1 



(41) 



The tilde denotes the imposition of the Dirichlet boundary condition at the outset. This 
expression yields an attractive force between the two Dirichlet points, 



F(a) 



dEo 



d(2a) 



dt 



i 1 



vr Vt 2 - m 2 (e 4at 



1) 



(42) 



In the massless limit it simplifies considerably: E 2 (a) = — n/48a and F(a) = — n/96a 2 . 
But this result is not internally consistent: as a — > oo, E 2 {a) ~^ 0, indicating that 
the energy of an isolated "Dirichlet point" is zero. The limit a — > also describes a 
single Dirichlet point, but E 2 (a) — > 00 as a — > 0. Also note that E 2 {a) is well defined as 
m — > 0, even though scalar field theories in one space dimension become infrared divergent 
when m — > 0. 

The Dirichlet point problem can nicely be studied with the presented method. A 
single delta-function background 



<j\ = XS(x 



has the Green's function 



G\{x,y) = G (x,y) 



XGo(x,a)G (a,y) 
1 + \Go(a, a) 



(43) 



(44) 



where Go(x,y) is the Green's function in the non-interacting case. The momentum argu- 
ment has been omitted. Obviously the limit A —>■ 00 gives Dirichlet boundary conditions, 
Goc^x, a) = Gqo (a, y) = 0. The vacuum polarization energy associated with eq. ([43]) is [H] 



£i(A) 



dt 

1,! 27T 



1 + 



21 



(45) 



Vt 2 - m 2 

To study the problem of two Dirichlet points it is obvious to consider 

a 2 (x) = X [5(x + a) + 8{x - a)} . (46) 
The renormalized Casimir energy for this potential has also been computed in Ref. ||, 

dt 1 



E 2 (a,X) 



2tt^ 



tin 



1 + 7 + 5 (1 " 



X 



(47) 



For any finite coupling A, the inconsistencies noted in E 2 (a) do not afflict E 2 (a, A): as 
a — > 00, E 2 (X) — > 2Ei(X), and as a — > 0, E 2 (a,X) — > Ei(2X). Also E 2 (a, A) diverges 
logarithmically in the limit m — > as it should. The force, obtained by differentiating 
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Figure 2: The energy density for the two delta function background at a = 0.2/m computed from 
eqs fll2"D and (27). Left panel: A = 3m. The distinct contributions associated with eqns. 
and (|27|) are disentangled but the singular contributions at x = ±a are omitted. Right panel: 
Dirichlet limit A — > oo in comparison to the the boundary condition result (|48|); A = (a — x)/a. 



eq. ( f47|) with respect to 2a, agrees with eq. ( (42|) in the limit A — > oo. However E 2 (a, A) 
diverges like A log A as A — > oo. Thus the renormalized Casimir energy diverges as the 
Dirichlet boundary condition is imposed, a physical effect which is missed in the pure 
boundary condition calculation. 

The counterterms to the energy density vanish away from x = ±a because they are 
local functions of a(x). Therefore the Casimir energy density for x ^ ±a can be calculated 
assuming Dirichlet boundary conditions from the start simply by subtracting the density 
in the absence of boundaries without encountering any further divergences 
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e 2 {x, a) 



e 2 (x,a 



oo 



m 

8a J in rr 

2 

m 



dtV^ 



m 2^cos T -f(x-a) 



4at 



iiu x ^ 

^ U n=l 



for \x\ < a 



H^)2 +m 2 



K (2m\x — a\) for |a;| > a . 

2ir 



(41 



This result excludes the points x = ±a. For finite A the Casimir energy density, e 2 (x, a, A), 
was also computed in Ref. || and is displayed in Fig. ^. The energy density between the 
isolated points is negative and approaches the boundary condition limit in a non- 
uniform manner. In the limit A — > oo it agrees with eq. ( pj| ) except at x — ±a where it 
contains extra, singular contributions. If one integrates eq. fl4"8| ) over all x, ignoring the 
singularities at x = ±a, one obtains eq. (|4lD . Including the singular contributions at ±a 
by integrating e 2 (x, a, A) gives eq. (|47|) . 

This simple example illustrates the principal situation: In the Dirichlet limit the 
renormalized Casimir energy diverges because the energy density on the "surface" , x = ±a 
diverges. However the Casimir force and the Casimir energy density for all x ^ ±a remain 
finite and equal to the results obtained by imposing the boundary conditions a priori, 
eqs. © and ©. 



Two Space Dimensions 

A scalar field in two space dimensions constrained to vanish on a circle of radius a presents 
a more complex problem. For smooth backgrounds only the local tadpole diagram diverges 
and thus the no-tadpole renormalization condition is still sufficient to render the theory 
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finite. For a(x) = X5(r — a) and r ^ a the subtracted local spectral density [pi(it,r)} 
vanishes exponentially as the momentum along the cut in eq. ( |T2"D increases. For finite A, 
both the t-integral and the £-sum are uniformly convergent so A — > oo can be taken under 
the sum and integral. The resulting energy density, e(r, A), agrees with e(r), obtained 
when the Dirichlet boundary condition, 0(a) = 0, is assumed from the start. As in one 
dimension, nothing can be said about the total energy because e(r) is not defined at 
r = a, but unlike the one dimensional case, the integral / dr e(r, A) now diverges even in 
the sharp limit for finite A. 

To understand the situation better, let's consider o~(x ) to be a narrow Gaufiian of width 
w centered at r = a and explore the sharp limit where w — > and o~(x) — > X5(r — a). 
For w 7^ 0, a does not vanish at any value of r, so [pe{it, r)] no longer falls exponentially 
at large t (and r ^ a), and subtraction of the first Born approximation to pi(it,r) is 
necessary, i.e. N = 1 in eq. (|l~2"l). As noted above, the compensating tadpole graph 
can be canceled against the counterterm, ciXa(x). The result is a renormalized Casimir 
energy density, e(r, w, A), and Casimir energy, E(w,X) = dr e(r, u>, A), both of which 
are finite. However as w — > both e(a, w, A) and E(w, A) diverge. 

The divergence can be traced to the 0(X 2 ) Feynman diagram. This diagram is sepa- 
rated by subtracting the second Born approximation to pe(it,r), i.e. N = 2 in eq. (|P2|). 
Then the £-sum and t-integral no longer diverge in the sharp limit, w — > but the equiv- 
alent diagram must be added back explicitly. In the limit w — >• it contributes 

A 2 a 2 r M 



/ dpJnfap) arctan (49) 

Jo 2m 



to the total energy. This diverges logarithmically as M — > oo. The divergence originates 
from the high momentum components in the Fourier transform of c(x) = X5(r — a) 
rather than the high energy behavior of the loop integral. This divergence gives an 
infinite contribution to the stress because it varies with the radius of the circle. This 



divergence only gets worse in higher dimensions (in contrast to the claim of Ref. |15|). 
For example, for cr(r) = X5(r — a) in three space dimensions the renormalized two point 
function is proportional to X 2 a 4 J M dpf(p) with f(p) = p 2 jl{pa)hip for large p. This 
integral diverges like Mln(M). 

The imposed upper limit, M, in eq. (|49| ) plays the role of a physical cutoff that regulates 
divergences localized on the surface. It is not related to the regulator of the ultraviolet 
divergences in loop integrals. Hence divergences like in eq. P9| ) are not renormalized by 
standard counterterms whose (divergent) coefficients are independent of the considered 
background because they are fixed by renormalization conditions on Green's functions at 
some prescribed finite external momenta. Divergences that emerge as M — > oo indicate 
that even the sharp limit w — > does not exist. The strong coupling limit A — > oo makes 
the divergence even worse. If the divergent terms depend on the quantity conjugate to the 
force (tension) under consider at ionP] the force (tension) cannot be defined independently 
of the structure of the material. This is the case for the shell or the sphere, but not for 
rigid bodies. 



3 For example, the distance between plates or the radius of the shell. 
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Conclusion 



In this talk I have presented an efficient method to compute vacuum polarization energies 
in renormalizable quantum field theories for static background fields. Starting point 
for this method is the energy density operator. Its matrix element in the ground state 
is expressed in terms of the Green's function which subsequently is parameterized by 
data from scattering off the background field. To compute the momentum integrals as 
contour integrals in the upper-half plane, this approach makes ample use of the identity 
of Feynman diagrams and Born approximants to Casimir energies. More importantly, 
this identity allows one to implement standard (perturbative) renormalization conditions 
on the divergent, low order Green's functions. In this way the removal of the ultraviolet 
loop divergences is independent of the considered background. 

Utilizing a variational approach to the so-obtained total energy, solitons can be con- 
structed. As an application I have shown that in a 1 + 1 dimensional chiral model quantum 
corrections create a soliton that is classically unstable. In similar 3 + 1 dimensional models 
the situation is more complex 0, as issues like Landau poles [IB] and sphaleron barri- 
ers [[H]] complicate matters. 

The divergences that arise when a quantum field is forced to vanish on a surface can 
be nicely studied in this approach by implementing a boundary condition as the limit 
of a less singular background. Energy densities away from the surfaces or quantities like 
the force between rigid bodies, for which the surfaces can be held fixed, are finite and 
independent of the material cutoffs. Observables that require deformation or change in 
area of the surface cannot be defined independently of the other material properties. 
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